Transients in Grounded Wires Lying on the 
Earth's Surface* 

By JOHN RIORDAN 

Voltages during transient conditions in a grounded wire lying on the 
earth's surface due to current in a second grounded wire also on the earth's 
surface are formulated for types of transient currents ordinarily obtained 
in a.-c. and d.-c. circuits. The fundamental formula is for voltage due to 
a unit step current, that is, a current zero for time less than zero, and unity 
for time greater than zero: curves are given for the function determining 
this voltage for a wide range of values of its two parameters. The for- 
mulas for other types of currents are not well adapted for numerical com- 
putation, which should be more conveniently carried out by numerical 
integration using the above curves. 

I 

A FORMULA for the mutual impedance of grounded wires lying 
on the earth's surface has recently been published by R. M. 
Foster. 1 The object of the present paper is to derive formulas for the 
voltages during transient conditions in one such grounded wire due to 
current in a second for types of transient currents ordinarily obtained 
in a.-c. and d.-c. circuits, and particularly for the voltage due to unit 
step current, zero for time less than zero, unity for time greater 
than zero. 

The voltage due to unit step current is expressed in closed form for 
straight parallel wires; closed form expressions have not been obtained 
for straight parallel wires for the exponential forms of current for 
a.-c. and d.-c. transients. While the integrals might be evaluated 
numerically, or transformed to asymptotic expressions, it appears 
more desirable in practical calculation to use the curves given for the 
unit step voltage directly; a single integration is necessary to find the 
voltage for current of arbitrary wave form, from the unit step result. 
The fundamental physical assumptions upon which the steady-state 
formula is based are as follows: The surface of the earth is assumed 
flat, the earth semi-infinite in extent, of uniform conductivity X, unit 

* A brief report of the results in this paper was given at the Summer Convention 
of the American Institute of Electrical Engineers, Toronto, Ontario, Canada, June 
23-27, 1930, in Discussion of "Mutual Impedances of Ground Return Circuits — 
Some Experimental Studies," by A. E. Bowen and C. L. Gilkeson; A. I. E. E. Trans., 
Oct. 1930. . 

1 R. M. Foster: "Mutual Impedances of Grounded Circuits" (Abstract), Bulletin 
of the American Mathematical Society, May, 1930, pp. 367-368; "Mutual Impedance 
of Grounded Wires Lying on the Surface of the Earth," Bell System Technical Journal, 
July, 1931. 
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permeability and negligible dielectric constant. The air above the 
earth is of zero conductivity, unit permeability, and negligible dielectric 
constant. Because of the assumption of negligible dielectric constant, 
the formulas for voltages during transient conditions do not hold 
strictly for small values of the time, that is, during the initial stages of 
the transient. The wires are of negligible diameter, lying on the 
surface of the earth, and insulated from it except at the ends, where 
there is point contact. 

In using the steady-state solution as the basis of transient solutions, 
the Heaviside operational calculus is employe d af ter replacing iu, 
where w = 2wf is the radian frequency and i = V — 1, by p = d/dt, the 
time differentiator, since (d n /dt")(exp icot) = (io)) n exp itat, where n is 
integral. 

II 

The mutual impedance of grounded wires lying on the surface of 
the earth and insulated from it except at the ends is given by the 
following formula: 2 

The integration is extended over the two wires 5 and s, having 
arbitrary paths, r and e are the distance and angle, respectively, be- 
tween differential elements dS and ds, and 7 = (4ttXuo) i/2 ; X is the 
ground conductivity and u = 2irf is the radian frequency. 

Replacing iu by p = d/dt in 7, the resulting forms to b e evaluated 
are exp ( — aV/>) and V/> exp ( — aVp) where a = rV^TrX. The first 
of these is known and, following Heaviside, 3 may be developed as 
follows. 

Expressing the exponential in series form: 

r--, rr , » 2 P a 3 p4p . 
exp(- «V/>) = 1 - aV£ 4--^ j! + '"• 

Integral powers of p are neglected, since (omitting the discontinuity 
at / = 0) the operand is unity and the derivative of a constant is 
zero. Then: 

« P (-«V?)-i-«V?[i+^ + #+•••]■ 

The bracketed terms may now be assumed to operate on Vp = (ttI)' 112 

2 Foster, loc. cit. 

3 Heaviside: "Electromagnetic Theory," Vol. II, pp. 49-51, equations (4) and (12). 



422 BELL SYSTEM TECHNICAL JOURNAL 

and, if p n is replaced by d n /dt", 

L 1 3*1! VW^ 5*2!^ JV^ 

since the term in brackets with its accompanying multiplier is the 
absolutely convergent expansion of the error function (erf) ; 



2 C* 

erf (z) = —= j exp ( — z 2 )dz. 

VirJo 



The result may also be established either by use of an integral 
equation 4 or the Fourier integral; it is given as pair 803, Table I, in 
tables published by G. A. Campbell. 6 In the present use of the tables, 
for unit step current, the mate of F(p)fp, where F(p) is a function of p 
to be evaluated, is taken since the unit step function is expressed by 
p- 1 (pair 415). 

The second operational form required may be derived from the 
first by differentiating with respect to a, since (d/da)F(p) = (d/da)f(t) 
where F(p) and /(/) are corresponding functions of p and /. Thus, 

a4p exp (- ajp) = ~ exp ( - ~ j , 
since 

| erf [>(/)] = JU'W exp | - [*(*)]* } ■ 

The unit step voltage may now be expressed, by substitution of 
these results, by the following formula: 

In equation (1), as in the steady-state formula from which it is 
derived, the wires are unrestricted in path or length on the surface of 

4 J. R. Carson: "Electric Circuit Theory and The Operational Calculus," McGraw 
Hill Co., 1926, p. 19, eq. 29. 

'"The Practical Application of the Fourier Integral," Bell System Technical 
Journal, October, 1928. 
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the earth. The formula for straight parallel wires, wire 5 extending 
along the z axis from — a to + a, and wire 5 from Zi to z 2 at distance 
x from it, is obtained by double integration between these limits with 
r 2 = x 2 + (5 - s)\ cos e = 1. 

The result of integrating once, with respect to S, is: 

7 "^ = 2^xJ 1 Ts L V* 2 + (a - s) 2 ~ V^ + (a + *) 2 J 

+ <t>(s + a) - 0(.v - a) J <fc, (2) 
where 

««) = . 5 erf ( V^T* J^ ) 

-^ exp (" Z T :! ) erf ("^)' 

where « is to be replaced by s -\- a and s — am equation (2). 

Equation (2) is checked as follows. In the first term substitute 
limits after removing differentiation and integration with respect to 
S, which cancel each other. In the second term integrate by parts: 



f \ erfj— IV+ (S - sYldS 

J r.v 2 + (s - s) 2 r \« L tw J J 



J [*« + {S - s) 2 _\ _ _ 

- , ^ c !/-»r l-7y^-«) a 

The integral coming from this operation combines with the remaining 
term to give : 

- *>f /?«*{- t [ * + {s ~ s)2 A dS ' 

which can be simplified in terms of the error function to the form in 
equation (2). 

Integration from Zi to z 2 gives the result: 

Vu(t) - tVdK** + °) - *(** - °) - ^ Si + °) + ^ Sl ~ a ^> < 3) 

where 

,/ \ * i V.r 2 + u- c ( r-n—. — 2 r^\ 



11 I 7rX.r 2 \ c I lir\ 

* ex p(--r) erf l"\T 
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As before, u is to be replaced in the equation by the functional 
arguments, which are the four sums of the s-coordinates of position. 
The factor x in \f/(n) is introduced to make it a function of two 
parameters, uxr 1 and ttXx 2 / -1 ; the result of integration is x~ l \f/(u). 
The result has the dimensions of abohms when all quantities are in 
electromagnetic c.g.s. units. 

To check equation (3) notice that the integration of the first term 
of equation (2) is effected by removal of differentiation and integra- 
tion signs, and substitution of limits; its contribution is identical with 
the d.-c. mutual resistance. 6 The integration of <£(«) may be effected 
by integrating the first term by parts and employing the indefinite 
integral: 



/ 



erf (ax)dx = x erf (ax) -\ p exp ( — a 2 x 2 ) + const. 

a-\T 



The result is checked by differentiating, that is, by the relation : 
%- \xr l +(u) + ; 1 1 = *(«). 

dill V.T 2 + U 2 J 

For large values of u, 

1 1 \ i u i r 1 / ^x 2 \ l 

*M~ — [l-exp^-— j , 

since 

erf (± oo) = ± 1, 

so that for a = °o the unit step voltage approaches the limit : 

-i[— (-mi 

where / = z 2 — Z\ is the length of the second wire. 

This result is in agreement with a result published by F. Ollendorff, 
Elektrische Nachrichten — Technik, October, 1930, eq. (26), and by L. C. 
Peterson, Bell System Technical Journal, October, 1930, equation (5). 

The case of collinear straight wires is obtained by taking the limit 
x = 0, which gives 



lim x V( w ) = - 
x=Q u 



- 1 + U + ~r rf rVT 



. /ttX / irX« 2 \l 

+ w VT exp (~~r/J 

= ur l ${u). 

This result involves the evaluation of an indeterminate form. 

9 G. A. Campbell: "Mutual Impedances of Grounded Circuits," Bell System 
Technical Journal, October, 1923, eq. (3), p. 5. 
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Fig. 1 — \p(u) for the range in which \p{u) < 1, ^ ux l ^ 10. 
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40 50 60 

VALUES OF UX"' 

Fig. 2 — >j/{u) for the range in which \p(ti) < 1, 10 < ux~ l < 100. 

Curves for \p(u) as a function of ux~ l with t/(Tr\x 2 ) as parameter of 
the curve families are shown on Figures 1, 2, and 3. The range 
\f/(u) — 1, is shown on Figures 1 and 2 for ux~ ] — 10 and 100, respec- 
tively; both figures cover the entire range of t/iirXx 2 ) in the intervals. 
The remaining range yp(u) > 1 is shown on Figure 3. For the greater 
part of the range on Figure 3 the function is determined by its limiting 
form for ux~ l large, that is, by the equation 



\f/(u) = ux~ l 1 — exp ( j— J 
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log iff(u) = log ux~ l + log 1 - exp ( - ^j- j • 

Thus Figure 3 may be used to indicate the range of applicability of 
the limiting form, which is quite large; in this range the unit step 
voltage is simplified as shown above. 
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Fig. 4 — The function f («). for collinear straight wires; for values below the range 

, . 1 , ttX?/ 2 
shown f («) ~ - - H — ■ 

The function £(«)> for the case of collinear straight wires, is shown 
on Fig. 4 for values of the argument //(ttXm 2 ) from 0.1 to 1000; for 
small values of the argument, the function is approximately 



. v 1 , ir\ll 2 



rXM 2 



s < 0.4 



These curves may be employed to obtain voltages due to other 
forms of disturbing currents by numerical or mechanical integration 
of the following integral: 7 



EM =J f P I(r)Vi 2 {t - r)dr 

= j f I{t - r)V 12 (r)dr, 
Jo 



where /(/) is the disturbing current as a function of time. 
7 J- R. Carson: loc. cit., p. 16, eq. (20) and (20a). 
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III 

The equation above may be used to obtain a formula for voltage 
due to suddenly applied current exp iut; or the operational product, 
of which it is an expression in terms of /, may be carried out directly 
in terms of p. The current is expressed in terms of p by: 

P 
exp iut = 



p — iu 



The second term in \j/(ii) is transformed by the operational equivalent 
already developed: 

erf — = = 1 — exp (— aV7>). 
2\7 

The last term in ^(«) is not known in closed form in p. 

The operational product of exp iut and the second term is evaluated 
by 

p[l - exp(- «Vft)l = P _ P exp ( - a-Jp) 
p — iu p — iu p — iu 

= exp iut — -z exp (/to/ — a\7to) erfc I — p — Vuo7 J 

rr- I a. rr— 

+ exp {iut + aV/co) erfc ( — p + V*'w/ 

the last term of which is given by pair 819 (with = 0) in the tables 
referred to. Erfc is the error function complement; 

erfc (s) = 1 — erf (2). 

The operational product of exp iut and the last term in \f/(u) may 
be expressed in integral form by the formula : 

-^-/W = [1+^1/(0 

p — 1U L p — lu } 

= f(t) + iu exp iut I exp (— iut)f(t)dt. 

The complete expression for the voltage due to cisoidal current is 
as follows: 

En(t) = tV [*( s 2 + fl ) - *( S2 - °) - *(* + °) + *<>< - «01 ( 4 ) 
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where 

u 2 u I t\x 2 \ ( jwK \ 

- "^'^ - [exp (iwt - 7 V.v 2 + « 2 ) erfc ( Jy (.t 2 + w 2 ) - JudJ 

+ exp (*w/ -f yV* J + w 2 ) erfc ( -*/— (.v 2 + m 2 ) + V7co7 I 

« . . , r* ( ■ tX* 2 \ .1 /irX \ j 
-twexpto)/ I exp ( — tut — I en I 11 \l— I at. 

The integral appearing in $(«) apparently cannot be expressed in 
closed form in terms of known functions; for numerical results series 
or asymptotic expressions may be derived but it appears more desirable 
to employ numerical or mechanical integration using the unit step 
voltage since tables or charts of the error function of complex variable 
which also appears in $(«) are not available. 

A useful check on the above formula is obtained by taking the limit 
for / = co , which gives the steady-state mutual impedance between 
straight parallel wires; the result is as follows: 

Z12 = £12(0 exp (— i<at) 
1 



27rX.r 
where 



0(s 2 + a) - <P(z 2 - a) - ¥(2, + a) + *(*i - a)], (5) 



V.x- 2 + ir 



ty(u) = — , exp (— yyjx 2 -+- u 2 ) 

xyx 2 -j- u 2 x 



Vx 2 + ir 



U f* ( y*X*\ . yu 

I exp — w . — erf — == dw 

* X \ 4w / 2Jw 

V* 2 + a 2 r . , j 2 , 2 , "I 

1 — exp ( — 7 Vx 2 + ii 2 ) 

iu r u , . 

- — I exp ( - 7 V* 2 + v?)dw. 



where as before 7 2 = Avkua. 

The t hird term in <£ («) approaches the limit given because 
erfc (— V« °°) = 2, erfc (V* °°) = 0; the integral term as given in 
the first form of ^(u) has been transformed by the substitution 

w = tut. 
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The first form of \P(m) may be checked directly from equation (3) 
by introducing iw = p in the operationally equivalent function of p; 
the third term of (3) being expressed by the infinite integral : 



F(p) =p f 'e-#f(t)dt. 
Jo 



The second form of ^(u) is obtained by separating the d.-c. mutual 
resistance term, and transforming the infinite integral as follows: 
express the error function in integral form, put y = yvj{2-4w) where 
y is the variable of integration for the error function, and invert the 
order of integration; thus 



exp — w — ei f — 1= aw 

V?r Jo Jo 



exp I — w — 



k 2 4- v 2 ) \ dw 
4w / Vw 



"V71 Jo Jo 



4s 2 



= 7 I exp ( — 7 V.y 2 4- i> 2 )Jy. 
Jo 



The infinite integral evaluated in the third line is No. 495 in Peirce's 
" Short Table of Integrals," third edition. 

The second form of >!>(«) may be verified by direct double integra- 
tion of the mutual impedance; it agrees with the known result in the 
limit for one wire infinite, and, when expanded in powers of 7, with 
the terms given in the second form for the mutual impedance by R. M. 
Foster, loc. cit. 

Expressions for voltages due to suddenly applied currents 
exp (— kt) sin to/ or 1 — exp (— kt), which are important forms for 
a.-c. and d.-c. networks, may be readily obtained from equation (4), 
the first by use of the expression: 

exp (— kt) sin ut = y= [exp (— kt + iut) 4- exp (— kt — iut)~\ 

and the second by the substitution — k = i<a and subtraction from 
the unit step voltage. 

The results attained in this paper depend in appreciable measure 
on advice and suggestions received from Mr. R. M. Foster of the 
American Telephone and Telegraph Company; I am also appreciative 
of the interest and advice of Messrs. K. L. Maurer and H. M. True- 
blood of this company. 



